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Effects of trace elements on the crystal field parameters of Nd ions at the 
surface of Nd2Fei4B grains 

Yuta Toga/ Tsuneaki Suzuki/ and Akimasa Sakuma^’'^^ 

Department of Applied Physics, Tohoku University, Sendai 980-8579, Japan 

Using first-principles calculations, we investigate the positional dependence of trace elements such as O and 
Cu on the crystal field parameter A2, proportional to the magnetic anisotropy constant of Nd ions placed 
at the surface of Nd2Fei4B grains. The results suggest the possibility that the parameter of Nd ions at 
the ( 001 ) surface of Nd2Fei4B grains exhibits a negative value when the O or Cu atom is located near the 
surface, closer than its equilibrium position. At the ( 110 ) surface, however, O atoms located at the equilibrium 
position provide a negative A2, while for Cu additions A2 remains positive regardless of Cu’s position. Thus, 
Cu atoms are expected to maintain a positive local of surface Nd ions more frequently than O atoms 
when they approach the grain surfaces in the Nd-Fe-B grains. 


I. INTRODUCTION 

Rapidly increasing demand for efficient electric motors 
motivates the development of high-performance Nd-Fe-B 
magnets. In sintered Nd-Fe-B magnets, Nd is frequently 
substituted with Dy, owing to Dy’s suppression of coer- 
civity (He) degradation. However, Dy is expensive and 
decreases the magnetization of Nd-Fe-B magnets. To 
realize Dy-free high-performance Nd-Fe-B magnets, we 
must understand the H(, mechanism of rare-earth (Re) 
permanent magnets. Although the He mechanism has 
been discussed for sintered Nd-Fe-B magnets,our un¬ 
derstanding of He is incomplete, requiring further exam¬ 
ination. 

Recent papers emphasized that the development of 
high-coercivity Nd-Fe-B magnets requires understanding 
their microstructure, especially the grain boundary (GB) 
phase surrounding the Nd2Fei4B grains.^’® Researchers 
have reported that the intergranular Nd-rich phase in¬ 
cludes neodymium oxides (NdO^,) with diverse crystal 
structures (e.g., fee, hep), suggesting that O atoms must 
exist near Nd atoms at the interface between the Nd-rich 
phase and Nd2Fei4B grains. Recently, Amin et al.^ 
confirmed the segregation of Cu to the NdOa;/Nd2Fei4B 
interface after annealing, suggesting that the magnetic 
anisotropy constant Ku grains are lapped by a Cu-rich 
layer. Thus, O and Cu atoms are considered to stay 
around the Nd ions at the interfaces, playing a key role 
for the coercive force of Nd-Fe-B sintered magnets. 

From theoretical viewpoint, it is widely accepted that 
the magnetic anisotropy of rare earth magnets is mainly 
dominated by the 4 f electrons in the rare earth ions, be¬ 
cause of their strongly anisotropic distribution due to 
their strong intra-atomic interactions such as electron- 
correlation and L-S coupling. The one-electron treat¬ 
ments based on the local density functional approxima¬ 
tion have not yet been successfully reproducing this fea¬ 
ture in a quantitative level, in contrast to the 3 d elec¬ 
tronic systems. To overcome this problem, crystalline 
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electric field theory combined with the atomic many- 
body theory has been adopted by many workers to study 
the magnetic anisotropy. In 1988 , Yamada et al.^ success¬ 
fully reproduced magnetization curves reflecting mag¬ 
netic anisotropy of the series of Re2Fei4B, using the crys¬ 
tal field parameters A™ as adjustable parameters. In 
1992 , the first principles calculations to obtain the A)" 
were performed by Richter et al.^^ for ReCos system and 
Fahnle et for Re2Fei4B system. 

Based on this concept, Moriya et al. showed via first- 
principles calculations that the crystal field parameter A2 
on the Nd ion exhibits a negative value when the Nd ion 
in the ( 001 ) planes is exposed to vacuum. As shown by 
Yamada et al.^ , the magnetic anisotropy constant 
originated from rare earth ion is approximately propor¬ 
tional to A2 when the exchange field acting on the 4 f 
electrons in a rare earth ion is sufficiently strong. Since 
the proportional coefficient is positive for Nd ion, nega¬ 
tive A2 implies the planar magnetic anisotropy. Further¬ 
more, Mitsumata et al} ’’ demonstrated that a single sur¬ 
face atomic layer with a negative value dramatically 
decreases He- In physical sintered magnets, however, the 
grain surfaces are not exposed to vacuum but instead 
face GB phases. Therefore, our next step is studying 
how interface elements in the GB phases adjacent to Nd 
ions affect the A2 values at the Nd-Fe-B grain surface. 
However, note that, in an actual system, many atoms in¬ 
teract with surface Nd ions and many possible configura¬ 
tions exist near the interface between GB and Nd2Fei4B 
phases. In this case, it is important from a theoretical 
standpoint to provide separate information on the posi¬ 
tional (r, 6 ) dependence of individual atoms on the local 
Ku (i.e., A2) of Nd ions at the grain surface. These anal¬ 
yses may provide useful information to judge the factors 
dominating the coercive force during experimental ob¬ 
servation of atomic configurations in real systems in the 
future. 

In this study, we investigate the influence of O and Gu 
atoms on the A2 of the Nd ion placed at the surface of 
Nd2Fei4B grains via first-principles calculations. As the 
sign of the single-site Ku of an Nd ion is the same as 
that of A2, the evaluation of A® is useful to understand 
how trace elements affect Ku at the surface of Nd2Fei4B. 
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FIG. 1. The geometric relationship between the Nd ion on the 
surface of Nd 2 Fei 4 B and the trace element for the (a) (001) 
and (b) (110) surface slab models. Here, r indicates the dis¬ 
tance between the Nd ion and the trace element, and 9 in¬ 
dicates the angle between the c-axis of Nd 2 Fei 4 B and the 
direction of r. This figure is plotted by using VESTA. 

We select O and Cu as trace elements on the basis of 
experimental results. 

II. COMPUTATIONAL DETAILS 

Electronic structure calculations were performed by 
density functional theory using the Vienna ab initio sim¬ 
ulation package (VASP 4.6). The 4f electrons in the 
Nd ions were treated as core electrons in the electronic 
structure calculations for the valence electrons, based 
on the concept mentioned in the previous section. The 
ionic potentials are described by projection-augmented- 
wave (PAW) method^'^ and the exchange-correlation en¬ 
ergy of the electrons is described within a generalized 
gradient approximation (GGA). We used the exchange- 
correlation function determined by Geperly and Alder 
and parametrized by Perdew and Zunger.* ’ 

We examined the (001) and (110) surfaces of Nd 2 Fei 4 B 
with the addition of the trace element using slab models. 
Figure 1 shows the geometric relationship between the 
Nd ion at the surface of Nd 2 Fei 4 B and the trace element 
for the (001) and (110) surface models. We placed the 
trace element at various distances r around the Nd ion 
at the surface, and at various angles 0 between the c-axis 
of Nd 2 Fei 4 B and the direction of r. In the (001) surface 
model (Fig. 1(a)), this unit cell has a vacuum layer equiv¬ 
alent to the thickness of the Nd 2 Fei 4 B unit cell along the 
c-axis (12.19 A ) and consists of 12 Nd, 58 Fe, and six B 
atoms. In the (110) surface model (Fig. 1(b)), we restruc¬ 
tured the Nd 2 Fei 4 B unit cell to expose Nd atoms on the 
(110) surface. The (110) direction in the original unit 


cell corresponds to the a-axis in the restructured unit 
cell. The restructured unit cell consists of 12 Nd, 68 Fe, 
and six B atoms. The lattice constant of the a-axis par¬ 
allel to the (110) surface is y/2 and perpendicular to it is 
l/y/2 times that of the original Nd 2 Fei 4 B unit cell. 

The (110) surface model has a vacuum layer with the 
thickness of 8 . 8/-\/2 = 6.22 A along the ( 110 ) direction. 
The mesh of the numerical integration was provided by 
a discrete Monkhorst-Pack fc-point sampling. 

To investigate the magnetic anisotropy of this system, 
we calculated A® for Nd ion adjacent to a trace element 
at the Nd 2 Fei 4 B surface. The value of the magnetic 
anisotropy constant originated from rare earth ion 
is approximately given by = —3J(J — l/2)a(r^)A2 
when the exchange field acting on the 4f electrons in a 
rare earth ion is sufficiently strong. Here, J is the an¬ 
gular momentum of 4f electronic system, a means the 
Stevens factor characterizing the rare earth ion, and (r^) 
is the spatial average of , as given by Eq. (3) below. 
For Nd ion, J = 9/2 and a is negative, and then posi¬ 
tive A 2 leads to positive Ku- The physical role of A 2 is to 
reflect the electric field from surrounding charge distribu¬ 
tion acting on the 4f electrons whose spatial distribution 
differs from spherical one due to the strong L-S coupling. 
Therefore, to obtain the value of A 2 , one needs the charge 
distribution surrounding the 4f electronic system with the 
following equationd'^di 

4 = -^^/ dRp[R)Zl{R) 


f 

^ / ^’’’’^^47rp4/(r)/(r2), 

(1) 

Z°2{R) = a{3Rl - \Rf)/\R\^ 

(2) 

(r^) = J (irr^47rp4/(r)r^. 

(3) 

where, a = (l/4)(5/7r)^/^, r< = min(r, i? ), 

and r> = 


max(r, |i2|). Here, P4f{r) is the radial part of the 4f 
electron probability density in the Nd ion, and p{R) is the 
charge density including the nuclei and electrons. The 
integral range of R in Eq. (1) is within a sphere with a 
radius of 70A from the surface Nd site. Here, the valence 
electron density was calculated by VASP 4.6, and P4f{r) 
was calculated using an isolated Nd atom. In p{R), the 
core electrons (including the 4f electron) forming pseudo¬ 
potentials in VASP were treated as point charges as well 
as nuclei. 


III. RESULTS AND DISCUSSION 

As shown in a previous paper, our computational 
procedure produces A 2 values of bulk Nd 2 Fei 4 B rea¬ 
sonably consistent with those obtained using the full- 
potential linear-muffin-tin-orbital method. In addition, 
our calculated A 2 for the Nd ion at the (001) surface of 
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FIG. 2. (a) The r dependence of crystal field parameter A 2 
and the total energy for O addition on the (001) surface of 
Nd 2 Fei 4 B for 0 = 0°. The closed squares and triangles indi¬ 
cate A 2 and total energies, respectively, (b) The 0 dependence 
of for r = 2.2 A, 2.5 A, 2.7 A, and 3.0 A. 


FIG. 3. Galculated valence electron density distributions for 
O addition on the (001) surface for r = (a) 4.0 A, (b)2.4A, 
and (c) 1.6 A, together with schematics corresponding to these 
three cases. 


Nd 2 Fei 4 B exhibits a negative value around —300K/a^ 
where ub represents the Bohr radius. This corresponds 
to that calculated using the full-potential linearized aug¬ 
mented plane wave plus local orbitals (APW-l-lo) method 
implemented in the WIEN2k code.^^ Therefore, our 
numerical calculation methods described in section II 
should be sufficiently reliable to quantitatively evaluate 
the influence of a trace element on acting on the Nd 
ion at the surface of Nd 2 Fei 4 B. 

Figure 2(a) shows variations in A 2 and electronic total 
energy when the O atom approaches the (001) surface 
Nd ion with angle 0 = 0° (see Fig. 1(a)). We confirm, 
as predicted by the previous work,^^ that A 2 exhibits a 
negative value when O is positioned at r > 3.5 A from the 
surface. At this situation, the total energy is confirmed 
to be almost equal to the summation of separated sys¬ 
tem of (OOl)-slab and O atom, each of which is -561.838 
and -1.826 eV. Thus the deviation energies from the sum¬ 
mation of these values represent the interaction energies 
between the slab and O atom. When O nears the Nd ion 
at 0 = 0°, A 2 becomes positive and increases to a peak 
value around 1000 K/a^ at r = 2.4 A. Interestingly, with 
further decrease of r, A 2 abruptly drops into negative 
values. 

To understand these behaviors, we show in Fig. 3 the 
calculated electron density distributions for r = (a) 4.0 A, 


(b)2.4A, and (c) 1.8 A, together with schematics corre¬ 
sponding to these three cases. The increase of A 2 with 
decreasing r for r > 2.4 A can be explained by the change 
in 5d electron distribution surrounding the 4f electrons 
of the Nd ion. That is, the 5d electron cloud extends 
towards the O atom through hybridization (Fig. 3(b)), 
which repositions the 4f electron cloud within the c-plane 
in order to avoid the repulsive force from the 5d electrons; 
this produces positive values of The decreasing be¬ 
havior for r < 2.2 A is attributed to the influence of the 
positively charged nucleus in the O atom exceeding the 
hybridization effect of the valence electrons (Fig. 3(c)); 
the attractive Coulomb force from O nucleus could rotate 
the 4f electron cloud so as to minimize the electro-static 
energy, resulting in negative values of A^. 

The variation of electronic total energy in Fig. 2(a) in¬ 
dicates a stable distance around r = 2.0 A, at which the 
value of A 2 is still positive. However, the value of A 2 eas¬ 
ily becomes negative with only a 0.2 A decrease from the 
equilibrium position, with an energy cost less than 1 eV. 
This deviation can take place due to stresses, defects, or 
deformations around grain boundaries in an actual sys¬ 
tem. This implies that the value of A 2 may exhibit a neg¬ 
ative value at real grain surfaces adjacent to GB phases, 
rather than to vacuum space. 

In Fig. 2(b), we show the 0-dependences of A 2 for 
various values of r. Starting from positive values at 
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O addition on the (110) surface 
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FIG. 4. (a) The r dependence of crystal field parameter A® 
and the total energy for O addition on the (110) surface of 
Nd 2 Fei 4 B for 6 = 90°. The closed squares and triangles indi¬ 
cate A 2 and total energies, respectively, (b) 0 dependence of 
for r = 2.2 A, 2.5 A, 2.7 A, and 3.0 A. 


FIG. 5. (a) The r dependence of crystal field parameter A 2 
and the total energy for Cu addition on the (001) surface of 
Nd 2 Fei 4 B for 0 = 0°. The closed squares and triangles indi¬ 
cate A 2 and total energies, respectively, (b) The 6 dependence 
of Ag for r = 2.7 A, 3.0 A, 3.5 A. 


0 = 0°, A 2 decreases monotonically with increasing 6 and 
reaches negative values for 6 > 50°. It is meaningful to 
compare these behaviors with those of the point charge 
model. In the point charge model, A 2 is proportional 
to Z 2 oc eq{3cos9^ — 1) (e > 0) where q is the valence 
number of the point charge. Thus the gross features of 
these results can be realized by the point charge model, 
if one assumes a negative charge {q < 0) for the O ion. 
Note, however, that we do not clearly identify the neg¬ 
ative charge on the O atom when r > 2.2 K. Therefore, 
the oxygen atom should be considered to influence the 
redistribution of valence electrons within the Nd atomic 
sphere such that the point charge model is applicable as 
if O were negatively charged. This can be regarded as a 
sort of screening effect.^® Similar explanations were pro¬ 
posed for the N effects in Re 2 Fei 7 N 3 ^'^ and ReFei 2 N,^^’^^ 
systems, where the N atoms changed the magnetocrys¬ 
talline anisotropy energies of these systems. 

Figure 4(a) shows the r-dependence of the values of A 2 
and electronic total energy when the O atom approaches 
the Nd ion in the (110) surface with an angle 6 = 90° 
from the c-axis. Contrary to the case for the (001) sur¬ 
face, A 2 is positive when Nd ions in the (110) surface are 
exposed to vacuum, and when the O atom is far from the 
surface. In this case, since the thickness of the vacuum 
space is only 6.22 A, the total energy exhibits symmet¬ 


ric variation with respect to r=3.lA. From this reason, 
we present total energies within the range r < 3.1 A. 
When the O atom approaches the Nd ion, the value of 
A 2 becomes negative for r < 3.5 A and stabilizes around 
r = 2.0 A, where the value remains negative. Notably, 
the value of A 2 becomes positive with only a 0.2 A de¬ 
crease in r. This is due to the attraction of the O nucleus 
to the 4f electron cloud, which aligns the 4f moment with 
the direction of the c-axis. The energy consumption for 
this reduction of r is less than 1 eV. In Fig. 4(b), we show 
the 0-dependence of the values of for varying r. As 
shown in Fig. 4(a), the values of A 2 for 6 = 90° are neg¬ 
ative in the range of 2.2 A< r <3.0 A. Similarly to the 
(001) surface in Fig. 2(b), A 2 increases with decreasing 
6. This can also be understood from the O atom’s redis¬ 
tribution of the valence electrons within the Nd atomic 
sphere, as if a negative ion exists in the direction of the O 
atom. The reason for the slight decrease in A 2 at 0 = 0° 
is not clear at this stage. 

Next we proceed to the case of Cu addition. Figure 5(a) 
shows the r-dependence of the values of A 2 and total en¬ 
ergy for Cu addition with 0 = 0°. The behaviors of A 2 
and total energy are almost the same as those for O addi¬ 
tion, shown in Fig. 2(a). However, the variations are less 
dramatic compared to the O case. This may reflect the 
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Cu addition on the (110) surface 
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FIG. 6. (a) The r dependence of crystal field parameter A® 
and the total energy for Cu addition on the (110) surface of 
Nd 2 Fei 4 B for 0 — 90°. The closed squares and triangles indi¬ 
cate A 2 and total energies, respectively, (b) The 6 dependence 
of for r = 2.7 A, 3.0 A, 3.5 A. 


weaker hybridization of the Nd atom with Cu than with 
O. The peak position of r (around 3.0 A) is greater than 
that with O addition, which may be due to the large 
atomic radius of Cu. In this case, a decrease of about 
0.5 A in r from the equilibrium position causes A 2 to be¬ 
come negative. This move costs around 0.3 eV in energy, 
which is less than that for O addition. The 0-dependence 
of A 2 for various values of r is shown in Fig. 5(b). The 
behavior can be explained through the hybridization be¬ 
tween valence electrons of both Nd and Cu atoms, as in 
the case for O addition. 

Figure 6(a) shows the r-dependence of the values of A 2 
and electronic total energy for the case of Cu addition to 
the (110) surface with the angle 9 = 90°. Contrary to 
the case of O addition in Fig. 4(a), A 2 maintains positive 
values for all r. This indicates weak interactions between 
Nd and Cu atoms. The steep increase of A 2 with decreas¬ 
ing r for r < 2.6 A reflects the significance of the nucleus 
effect of Cu on the 4f electron cloud. The weak hybridiza¬ 
tion between Nd and Cu atoms for r > 2.6 A can be seen 
in the 0-dependence of A^^ as shown in Fig. 6(b). 


IV. SUMMARY 

Motivated by a recent theoretical work'^ demonstrat¬ 
ing that a surface atomic layer with negative magnetic 
anisotropy constant can drastically decrease the co- 
ercivity He, we evaluated the influence of trace elements 
O or Cu on the crystal field parameter A 2 of the Nd 
ion, at both the (001) and (110) surface of the Nd 2 Fei 4 B 
grain. 

In both cases of O and Cu additions to the (001) sur¬ 
face, decreasing the distance r to the surface with con¬ 
stant 0 = 0° changes A 2 from negative to positive values. 
At the equilibrium position, the value of A® is found to 
remain positive. With further decrease of r, A 2 abruptly 
becomes negative again with an energy cost less than 
1 eV. This is due to the positive charge of the nucleus 
in O and Cu, which gains influence with decreased dis¬ 
tance from Nd. Therefore, the value of A 2 may exhibit 
a negative value due to stresses, defects, or deformations 
around grain boundaries adjacent to the GB phases in 
an actual system. 

The 0-dependence of A 2 can be roughly expressed as 
3 cos^ 0—1 for both O and Cu additions. This sug¬ 
gests that these elements redistribute the valence elec¬ 
trons within the Nd atomic sphere such that the negative 
point charge model is applicable, as if these species had a 
negative charge. We observed that the strength of these 
addition effects is larger with O than with Cu. This dif¬ 
ferent behavior of O and Cu atoms is clearly seen in the 
case of (110) surface. For 0-addition, the r-dependence 
of A 2 is opposite to that in the (001) surface case, as is 
expected from geometrical effects. Actually, the surface 
Ku potentially decreases to negative values at the equi¬ 
librium position of O. However, For Cu-addition in the 
case of (110) surface, the variation is small compared to 
O addition, and A® remains positive for all r. 

Therefore, O is expected to produce negative interfa¬ 
cial Ku more frequently than Cu when it approaches the 
Nd ion at the grain surface. The analysis of the total 
energy showed that local stable positions of the trace 
element exist for the special configurations considered 
here. However, due to the complex interatomic inter¬ 
actions and local stresses in real multi-grain structures 
of Nd-Fe-B magnets, many possible configurations exist 
in the local crystalline structure near the interface be¬ 
tween GB and Nd 2 Fei 4 B phases. In this sense, the (r, 
0) dependences of the local (i.e., A 2 ) shown in this 
study may apply when we consider the effect of individual 
atoms adjacent to Nd ions at the interfaces of GBs. 
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